function result = modelPred(XTRAIN,ytrain,XTEST)
b = glmfit(XTRAIN, ytrain,'normal');

%b = regress(ytrain,XTRAIN);
result =[ones(size(XTEST,1),1) XTEST]*b;


% size(ytrain)
% size(XTRAIN)
% size(XTEST)
% size(result)
% size(b)
%diff = ytrain' - result;
%tMSE = mean(sum(diff.*diff));
%  result = XTEST*regress(ytrain,XTRAIN);
% disp('test')
% 
% figure
% plot(t,result,tytrain,'LineWidth',1.8,'MarkerSize',8);
% legend('Estimiated Traj','Actual');
% title(['Errors: ' num2str(tMSE)]);

end